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ABSTRACT 

We observed the Termi-discovered 7 -ray binary IPGL J1018.6—5856 at 20 epochs over 50 days using the 
CHIRON spectrograph, obtaining spectra at R ~25,000 covering AA4090 — 8908A. The average spectrum 
conhrms an 06 V((f)) spectral type and extinction E{B — V) = 1.35 ± 0.04. Variable absorption line equiv¬ 
alent widths suggest substantial contamination by wind line features. The limited S/N ratio hindered accurate 
continuum dehnition and prevented measurement of a high quality radial velocity curve. Nevertheless, the 
best data indicate a radial velocity amplitude < 40km for the He II lines and substantially lower for H I. 
We argue that this indicates a most likely compact object mass < 2.2Mq. While black hole solutions are not 
excluded, a neutron star source of the 7 -ray emission seems preferred. 

Subject headings: gamma rays: stars — binaries: spectroscopic 


1. INTRODUCTION 

To date only hve binaries with high mass secondaries have 
been detected in the GeV band; IPGL 1018.6—5856 is the 
only one whose binary period was discovered at these ener¬ 
gies (ICorbet et al.ll20HT lAckerman et al.ll2012^ . The nature 
of the 7 -producing compact object in these binaries is un¬ 
clear: it might be a jet-producing black hole, in analogy to 
7 -loud blazars (the ‘fiQSO model’), or the 7 -rays may be 
generated by the magnetosphere or wind termination shock 
of an energetic pulsar (the ‘PSR/PWN model’). This latter 
model certainly applies to the PSR B1259—63 system, whose 
GeV emission is produced via Inverse Compton scattering 
off wind shock particles when the PWN is compressed near 
periastron. For IFGL 1018.6—5856 and the other three sys¬ 
tems (LS Ih- 61°303, LS 5039 and Cygnus X-3), the identity 
of the compact object is presently unknown. Discovery of 
7 -ray pulsations would certainly indicate a neutron star pri¬ 
mary. However, the poorly known acceleration in the binaries 
makes 7 -ray pulse searches prohibitively expensive and the 
dense plasma produced by the massive secondary winds frus¬ 
trates radio pulse searches. An alternative approach is to con¬ 
strain the compact object mass via measurement of the sec¬ 
ondary’s radial velocity. Although the high mass and strong 
winds of the secondary also make this very challenging, even 
crude constraints can be useful; for example a minimum com¬ 
panion mass > 3Mq would certainly indicate a black hole. 

Aft er discovery of a 16.5 d 7 -ray modulation (ICorbet et alJ 
I2OI 111 , the associated Swift X-ray source was identihed 
with a bright optical counterpart. A low resolution spec¬ 
trum with the SAAO 1.9 m shows an 06V((f)) companion 
with an estimated mass « 20 — 30Mq and distance « 5 kpc 
(lAckerman et al.ll2012^ . This system is also a radio source 
and Swift XRT (An et al 2013) and ATCA monitoring detect 
the orbital modulation in the X-ray and radio band, respec¬ 
tively. Optical photometric monitoring detects no variation, 
although present constraints are weak. These properties are 
similar to those of the more compact (Pf,=3.9 d) 7 -ray binary 
LS 5039. There spectroscopic and photometric studies (e.g. 
ISartv et ^1201 Ih indicate an 06.5V((f)) secondary, a sub¬ 
stantial eccentricity e = 0.24 ± 0.08, and a modest mass 
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function f{m) = 0.0049 ± O.OOOGMq. The limit on LS 
5039’s photometric variability, 2mmag, implies a high sys¬ 
tem inclination and/or a small compact object mass. Since 
the inclination is poorly determined, they can only conclude 
that Ml > 1.8Mq, allowing both black hole and neutron star 
solutions. 

In this paper we describe a project to study the radial ve¬ 
locity of IFGL 1018.6—5856 using echelle spec t roscop y. We 
adopt the orbital period estimate of IColev et al.l (12014!) Pb = 
16.549 ± 0.007d with epoch MJ D 55403.3 ± 0.4 corresp ond- 
ing to the 7 -ray peak described in lAckerman et al.l (l2012h . An 
independent Swift/Fenni analysis (An, priv. comm.) con¬ 
hrms this ephemeris. While our results appear to show sys¬ 
tematic variation with orbital phase, it is unclear if we have 
isolated the true radial velocity curve. Nevertheless, our mea¬ 
surements improve the characterization of the companion star, 
limit the amplitude of its motion and point to the sort of ob¬ 
servations needed for a dehnitive study. 

2. OBSERVATIONS AND CALIBRATIONS 

The 16.5 d binary period is particularly inconvenient, being 
too long to cover in a dedicated run, but too short for effec¬ 
tive coverage with ad hoc ^ monthly observations. The best 
solution was queue observing and we were awarded 25h of 
queue time (NOAO-1 lB-0092) on the CTIO 1. 5 m using the 
CHIR ON hber-coupled echelle spectrograph (ISchwab et alJ 
IMoh . IFGL 1018.6-5856 is relatively faint for CHIRON 
alV = 12.5, so CHIRON was run in the slitless (2.5" hber 
aperture) mode for an effective resolution R ~ 25, 000. The 
CCDs were binned 4 x 4 to maximize S/N and minimize read¬ 
out time. The spectra covered 4090—8908A, with substantial 
order overlap except in the reddest orders. 

Observations were attempted roughly every other night 
from March 02, 2012 (JD 2455989) to April 20, 2012 (JD 
2456022). Four nights were lost to system failures and 
weather; spectra were obtained on 20 nights (Table 1). We 
requested 3x 1000s exposure on IFGL 1018.6—5856, one 
night had only a single exposure. For each epoch we obtained 
ThAr arc exposures immediately before the target observa¬ 
tions for wavelength calibration. For the first epoch only we 
also observed the flux standard HD60753. CHIRON was un¬ 
dergoing final commissioning during these observations, and 
the obtained S/N was substantially lower than predicted by 
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TABLE 1 

CHIRON Observations 


Date 

Nexp 

S/N 

6000A 

4 >b 

03/02 

3 

36 

0.41 

03/06 

1 

24 

0.65 

03/10 

3 

12 

0.89 

03/14 

3 

33 

0.13 

03/16 

3 

32 

0.25 

03/18 

3 

44 

0.37 

03/20 

3 

37 

0.49 

03/22 

3 

33 

0.61 

03/24 

3 

3 

0.73 

03/26 

3 

31 

0.85 

03/28 

3 

18 

0.98 

03/30 

3 

18 

0.09 

04/02 

3 

28 

0.28 

04/04 

3 

31 

0.40 

04/06 

3 

21 

0.51 

04/08 

3 

8 

0.63 

04/10 

3 

22 

0.76 

04/12 

3 

20 

0.87 

04/18 

3 

30 

0.24 

04/20 

3 

29 

0.38 


the design plots (predicted peak S/N/extracted pixel ~ 40 in 
each 1000s exposure). In practice, we obtained S/N/extracted 
pixel ~ 25 at the order peaks in our best exposures (in line 
with recent revised 6.5% instrument efficiency). More typi¬ 
cal S/N/extracted pixel was Ri 17, with a few epochs much 
worse. We were not able to ascertain from the night logs why 
the S/N on source varied so dramatically, although from the 
low apparent flux at some epochs, we suspect that poor trans¬ 
parency and imperfect fiber acquisition may play a role. In the 
end we analyzed the 17 epochs with S/N> 15 in the summed, 
combined spectra. These spectra sample four binary orbits. 

Although the CHIRON data distribution included a stan¬ 
dard pipeline extraction of the echelle spectra, this was op¬ 
timized for the relatively red target stars of the instrument’s 
principal (radial velocity planet search) program. This ex¬ 
traction covered AA4504 — 8908A (62 orders); for our blue 
06V target 14 additional orders were measurable, extending 
the spectral coverage to 4090A. 

Standard CCD processing was performed with the IRAF 
package. A master bias was assembled and subtracted from 
the target, flat and lamp images. The spectra for each epoch 
cover 3200s. As this is only 0.2% of the orbital period, 
intra-epoch velocity smearing was negligible and the three 
target exposures for each night were median combined with 
average a-clipping to increase S/N and reject cosmic ray 
events. Combined flat field images were also prepared for 
each night. These were used to define the apertures for each 
order. Sufficient inter-order pixels were available to define 
flanking ‘background’ windows for each order. A linear 
fit to these backgrounds was used as a local scattered light 
subtraction. All extractions were performed using variance 
weighting, modeled by the CCD read noise (4.5e“) and gain 
(1.3e-/ADU). 

The extracted science spectra for each order were divided 
by the extracted flat spectrum. This removed the blaze func¬ 
tion and summed pixel-to-pixel variations, while maintaining 
a relatively stable response from the smooth lamp spectrum. 
In principle, the true response could be determined from the 
standard star spectrum. In practice, for the queue-selected 
standard the tabulated fluxes were derived from a spectral 
model and were too sparse (lOA grid) to fit an accurate sensi- 
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Fig. 1.— Fluxed CHIRON spectrum from March 02, 2012, corrected for 
atmospheric and interstellar (i?B—v = 1-35) extinction. Telluric features 
have not been removed. The strongest night sky and absorption lines are 
visible, along with inter-order gaps for the reddest five orders. 

tivity function for each order. Accordingly we used the lamp- 
flattened spectra directly. The wavelength calibration was de¬ 
termined from the nightly ThAr spectra. Non-linear solutions 
were fit, with a typical rms residual 0.01 A. Note that although 
CHIRON is designed to be quite stable, night-to-night shifts 
in the spectrum position on the chip were significant. Aper¬ 
ture migrations (1-2 pixels) were treated by re-centering or¬ 
ders on the flat before extraction (and importing these traces 
to the arc and science spectra). Shifts along the dispersion di¬ 
rection (~ O.lA^ 5 — 10kms“^) were removed by running 
IRAF task re-identify on each night’s ThAr arcs to derive the 
wavelength calibration. 

After trimming very low S/N data at the edges of each 
echelle order, we merged all orders to a single spectrum, 
average-combining the overlap regions (figure [l]i. Given our 
relatively modest S/N for most of the spectra, this proce¬ 
dure generally produced an adequate match between orders. 
For some epochs, small flux gradients across orders and poor 
inter-order match were however present. Such imperfections 
can be attributed to differences in illumination paths of and/or 
light path shifts between flat and science images. These may 
be removed by computation of a smoothly varying residual 
blaze function (iSkoda et al.l 120081) . Although we attempted 
such corrections, the S/N of our targets and standard expo¬ 
sures did not allow significant improvement. Finally, we di¬ 
vided the merged spectra by a low order function to remove 
the global lamp response and create normalized spectra ap¬ 
propriate to our kinematic/line study. The S/N weighted com¬ 
bined normalized spectrum from all epochs is shown in figure 
|2] An absolute fluxed spectrum was computed from the first 
night’s data where the adjacent standard was used to normal¬ 
ize the fluxes. For each epoch we corrected the normalized 
spectrum to the heliocenter using the IRAF routines rvcorrect 
and dopcor. 

3. SPECTRUM ANALYSIS 
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Fig. 2. — Normalized average spectrum. Each sub-panel covers lOOOA. 
The strongest lines used for the radial velocity analysis are labeled, addi¬ 
tional lines from these species are ma rked. The synthesized D IB transmission 
spectrum is generated from tables in IJenniskens et all <19961) (normalized to 
Eb-v = 1-35) and is plotted offset by -0.1. 

The overall spect rum (figure [B is quite c onsistent with the 
SAAO spectrum of I Ackerman^ al] (1201 2l) . with H I, He I 
and He II lines competing with strong interstellar absorption. 
The He IlA4541/He I A4471 ratio is a go od diagnostic of sub¬ 
class for O-type stars (iGrav et al.ll2009l) . From our combined 
spectrum we measure 2.3 ifc 0.2 for this ratio, indicating 06 
dJascheck & Jaschecldll995h . The N III emission line triplet 
4634/4640/4642 is not present, or even weakly in absorption. 
Even higher S/N is needed to measure the state of this line. 

Its weakness argues for a dwarf classification; overall the des¬ 
ignation 06 V((f)) is supported by this spectrum. The photo- 
spheric lines are broad, reaching FWHM 350 kms“^. There 
were no sky fibers for sky subtraction and so a few of the 
brightest night sky lines (e.g. 5577A ) are visible. A broad 
emission feature at 4293A appears in most of the spectra, with 
highly variable flux and high equivalent width even when the 
target was weak. As this corresponds to no strong astronom¬ 
ical line, we suspect a spectrograph artifact from local light. 

In figure|2]we label the strongest photospheric features. 

3.1. Interstellar Absorption 

Our measurement of the photospheric features is com¬ 
plicated by particularly strong interstellar absorption. 
Atomic features (e.g. Mg I 5174, Na I 5892) provide 
narrow absorptions, while Diffuse Interstellar Bands (DIBs) 
are particularly prominent in the spectrum . To identify 
these, we use the catalog of iJenniskens et all (119961 see also 
http://leonid.arc.nasa.gov/DIBcatalog.html) 
to synthesize a DIB absorption spectrum, combining Gaus¬ 
sian features at the cataloged Agg, kinematic and relative 
equivalent widths (figure|2]l. Comparison with the companion 
spectrum shows that many subtle broad features have a 
DIB origin. Note that the relative feature strength shown 
here is only typical; substantial variation exists between 
various Galactic sight-lines. The many DIB features make 


TABLE 2 

Reddening from DIB Measurements 


DIB 

A 


bt 

EW 

A 

Eb-v 

5780 

0.506 ± 0.050 

-0.014 ±0.032 

0.709 ± 0.009 

1.43 ±0.16 

5797 

0.124 ±0.014 

-0.003 ± 0.009 

0.134 ±0.006 

1.10 ±0.14 

6196 

0.048 ± 0.004 

0.001 ±0.003 

0.063 ± 0.005 

1.29 ±0.12 

6379 

0.057 ±0.012 

0.004 ± 0.008 

0.119 ±0.005 

2.00 ± 0.44 

6613 

0.192 ±0.026 

0.003 ±0.017 

0.225 ± 0.006 

1.16 ±0.18 

6660 

0.014 ±0.004 

0.009 ± 0.003 

0.043 ± 0.004 

2.54 ± 0.82 


^ Coefficients from IKos et al.l I2013D . 


it particularly challenging to define the true photospheric 
continuum, in turn this can subtly affect our line centroid 
estimates. 

The DIBs are believed to be generated by PAHs associated 
with interstellar dust. Thus one expects a correlation between 
DIB strength and interstellar reddening from dust scattering. 
iKos et al] (1201 3h lists fit linear relations for several strong 
DIBs of the form EW = a x Eb-v + b. These relations de¬ 
pend on whether the sight-line is ‘UV shielded’ (C sight-line) 
or not (a sight-line), and the A5797/A5780 EW ratio is used 
to discriminate these cases. Our value is 0.19 ± 0.01, below 
the 0.3 threshold, so we treat this as a cr sight-line. In tabled 
we apply the tabulated coefficients to several well-measured 
DIBs to estimate the reddening. The inferred Eb-v have a 
surprisingly large scatter, with a weighted average Eb-v = 
1.35±0.04 (bootstrap error). This value agrees v ery well with 
the est imate Eb-v = 1.34ib0.04 determined bv l^poli et alJ 
(1201 Ih from Stromgren photometry. It is similarly in reason¬ 
able accord with the Nh estimated from the X-ray spectrum. 
In turn, this agreement supports the photometric distance es¬ 
timate d = 5^.2 1 kpc. 

3.2. Radial Velocity Measurements 

To confirm the stability of our wavelength solutions, we 
measured the night sky lines A5577.3387 A and A6300.304 A 
in each processed spectrum, before heliocentric corrections 
were applied. The rms wavelength variations were 0.0106A 
and 0.0124 A respectively, i.e. aRv — 0.6kms“^. 

A second stability check was made measuring the very 
strong Na(D) interstellar doublet A5889.951/5895.924A in 
spectra corrected to the heliocenter. The line centroids dis¬ 
play mean redshifts of 2.88 and 2.85 kms“^, with spectrum- 
to-spectrum rms variation of 0.88 and 0.74kms“^, respec¬ 
tively. We conclude that our wavelength solutions are good to 
better than 1 kms“^. 

To attempt to measure the photospheric velocity of the 
companion, we must follow the centroids of the broad (^ 
300km s“^) H I, He I and He II lines, to a few kms“^. This 
is challenging since the imperfect inter-order fluxing causes 
significant gradients and steps in the continuum, with a char¬ 
acteristic scale comparable to the typical 50A (3000 kms“^) 
order width. Eurther, for several of the photospheric lines 
there are strong DIBs in the line wings, which can perturb 
the continuum. To mitigate this latter effect, we divide the 
spectrum by the model DIB spectrum. While imperfectly re¬ 
moving all DIBs, this does significantly flatten the continuum 
in the region of the strongest features. 

All photospheric lines were fit with simple Gaussian pro¬ 
files (test Voigt fitting showed negligible Lorentzian wings). 
The fitting used the IRAE deblending routine, where the ef- 
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Binary Phase (0^) 

Fig. 3.— Radial Velocity measurements for three species. Points with 
circles are from high S/N spectra; those with crosses from lower S/N data. 
Black squares show the mean velocity of the two strongest lines for H I and 
He II for the high S/N spectra; large systematic scatter is seen. The sample 
K = 40kms“^ radial velocity curve plotted provides a reasonable match 
to the He II average points; the H I amplitude is ~ 2.5X smaller (see text). 

fective noise was determined for each line by a measurement 
of the spectrum rms in the nearby, line-free continuum. The 
statistical errors on the fit parameters (line centroid, strength 
and, optionally kinematic width) are provided by the routine, 
based on Monte Carlo simulations. For the lower S/N lines the 
kinematic width was held fixed at the value measured from the 
average spectrum (5-lOA FWHM). 

We also attempted to determi ne the line centroids usin g 
the so-called “mirroring” method (iParimucha & Skodal2008^ . 
This allows one to separately probe the line wings and core. 
Although we can see that our lines are significantly asym¬ 
metric, suggesting absorption in wind/outflow components, 
we lacked the S/N to measure these structures accurately and 
compare them across lines and species. This prevented us 
from using the line shapes to improve our estimate of the sys¬ 
tematic (photospheric) radial velocity. 

4. RADIAL VELOCITY VARIATIONS 

To test for velocity variation, we performed a weighted 
least-squares fit for the mean velocity of each species. Table 
3 shows the velocities, errors, goodness of fit and the proba¬ 
bility of constancy. It is unclear how to compare mean veloci¬ 
ties between lines, as the effective rest wavelengths depend on 
the excitation. Also, the measured line centers are somewhat 
uncertain due to continuum irregularities; for example He II 
6527 appears substantially redshifted, due to a continuum gra¬ 
dient and superposed DIB. Finally, the lower excitation lines 
can be affected by w ind absorption and may be blue shift ed 
with respect to He II (ICasares et alJl2005l : ISartv et alJl201 Ih . 

Overall, we infer a net systemic radial velocity of ^ 
40kms“^, but a precise value requires cross-correlation with 
a simultaneously observed RV standard of similar spectral 
class. The reduced values are, however, interesting. From 
the low probabilities, we infer that the the strong lines do in¬ 
deed show statistically significant variation. 

However, there are certainly large systematic effects, and it 


TABLE 3 

Test for Radial Velocity Variations 


Species 

■^rest 

A in air 

r 

km 

X^/DoE 

Log(Prob) 

Ha 

6562.808 

43.0 ± 1.1 

14.7 

-3.9 

H/3 

4861.34 

22.4 ±3.0 

3.6 

-1.2 

H7 

4340.48 

25.5 ±36.7 

3.3 

-1.2 

He I 

5875.65 

72.9 ±4.3 

20.8 

-5.3 

He II 

4685.682 

44.1 ± 20.9 

10.7 

-3.0 

He II 

5411.52 

63.3 ±3.3 

15.5 

-4.1 

He II 

6527.099 

125.8 ±15.4 

5.6 

-1.7 
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Binary Phase (0j) 

Fig. 4.— Equivalent width and its variation for strong absorption lines. H/3 
and He 11 4686 absorption appears to weaken near phase 0^ = 0.2. 
is unclear how much of this variation is associated with the 
photosphere radial velocity. To illustrate, we average the ra¬ 
dial velocity from the two strongest HI lines (Her and H/3) and 
two strongest He II lines (A4686 and A5412) and plot the val¬ 
ues for the highest quality (S/Ngooo > 20) spectra (Figure 3). 
Even disegarding the mean radial velocity F offset between 
lines, we see that the fluctuations are several times the statis¬ 
tical error. Nevertheless, both sets of velocities show similar 
trends, with relatively low values in the range = 0.0 — 0.2. 
Since (pB = 0 is the 7 -/X-ray light curve maximum, it is 
tempting to identify this excursion with periastron and inter¬ 
pret the radial velocity measurements as showing significant 
eccentricity. The two species H I and He II show mean veloc¬ 
ities Fh I = 31.9kmand Fne ii = 46.2kmand RMS 
dispersions cth i = 11.1kmand dHe ii = 29.2kms“^. 

4.1. Equivalent Width Variations 

Following ISartv et ^ (1201 ll) . we checked the equivalent 
widths of our strongest absorption lines. Ha, the strongest 
atomic absorption has an average equivalent width of — 2.2A. 
Like these authors we also see evidence for equivalent width 
variation through the orbit, especially for H^ and He II 4686 
(Figure 4). However, in our case the minimum absorption 
strengths appear to lie at i/ib « 0 .2, close to compact object 
superior conjunction for the orbit suggested in figure [3] 0-e. 
when viewing the side of the companion away from the com- 
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pact object). This might indicate denser wind flow through 
the L3 point. Inspection of the strongest lines show variable 
wings, suggesting time or phase-dependent views of such out¬ 
flow. 

As for LS 5039, we find that our observed Ha equiv¬ 
alent width, —1.87 to — 2 . 29 A, is smaller than the ex- 
pected pa —3.27A ph otospheric width for this spectral class 
(iMcSwain et al.l[^04l) . Scaling from those authors’ estimates 
(i.e. assuming similar companion mass and asymptotic wind 
velocity), we find that our lower observed residual absorption 
equivalent width implies a somewhat larger mass flux, namely 
M « 0.7 - 1.0 X 10-^Mq y-i. 


4.2. Radial Velocity Model and Implied Mass 

Labeling the primary (compact object) 1 and the visible sec¬ 
ondary 2 , we have a single line spectroscopic binary, with ra¬ 
dial velocity measurements constraining the mass function 


/(Ml) = 


Mf 


(Ml -f M2Y 




with K 2 = AV 2/2 the radial velocity amplitude of the sec¬ 
ondary and i the unknown system inclination. Since sim < 1 
and M 2 > 0, f(M^) provides a lo wer bound to the mass 
of the primary. fCasares et al.l (l2005h give the mass of an 06 
V((f)) star as 20.0 - 26.4 Mq. 

In view of the large and poorly understood systematics, one 
should not directly fit our radial velocities. However, a ba¬ 
sic model (figure gj gives an idea of the allowed amplitude. 
Here a curve is plotted for K = 40km s“^, L = 40km s“^, 
e = 0.2, and longitude of ascending node w = O.Ott. This 
provides a reasonable match to the range of the He II line; the 
H I range is K Ki 15km s“^. Larger amplitudes or eccentrici¬ 
ties provide a poor match to the data, so we take th is as an up¬ 
per lim it on the allowed radial velocity amplitude. ISartv et ^ 
(1201 ll) adopt the He II radial velocity curve as minimizing 
wind contamination and best following the photospheric mo¬ 
tion. Their curve for H I has a somewhat smaller amplitude 
than that for He II; for our data t he range of H I ve locities is 
^ 2.5x smaller. Thus we follow l^rtv et al.l (1201 Ih in adopt¬ 
ing the He II velocities, but consider the possibility that these 
are systematics dominated and that the stronger H I lines re¬ 
flect the true radial velocity amplitude. 

Figure |5] shows the corresponding allowed regions in the 
mass-mass plane for circular orbits. For the allowed range 
of eccentricity Mi will be up to 7% smaller. The solid lines 
show the K = 40km constraints suggested by He II, the 
dashed lines show the range from the Balmer lines. 


5. DISCUSSION AND CONCLUSIONS 

Our attempt to measure the radial velocity curve of IFGL 
1018.6—5856 has been frustrated by a large scatter, appar¬ 
ently due to systematic errors. In part these are certainly due 
to the imprecision of the continuum calibration in our echelle 
spectra. In part they are also likely due to variable absorption 
line components of a stellar wind unpredictably shifting the 
line centroids. 

However, we do see that the maximum radial velocity am¬ 
plitude allowed by our data is ^ 40km s“^. If this is the true 
radial velocity amplitude, then we would require a compact 
object mass > 1 . 7 M 0 and 50% of the allowed parameter 
space has Mi > 2.2Mq. Since the likely maximum neutron 
star mass is 2 . 5 M 0 , this would seem to suggest that the 
7 -ray source is a black hole. However, one should compare 


with th e masses of other stellar mass black holes. iFarr et all 
(1201 ll) find that at 90% confidence the minimum Mbh in 
stellar mass binaries is > 4 . 5 M 0 . Although they do not sep¬ 
arately evaluate this for the high mass systems, like IFGL 
1018.6—5856, their figure 2 indicates most probable masses 
10 — 2 OM 0 and 90% lower mass limits > IOM 0 . From fig¬ 
ure |5] we see that we require i < 30° for Mi > A.^Mq {a 
priori probability 13%), while for Mi > IOM 0 we require 
i < 14° (probability 3%). Thus the likelihood of having a 
typical stellar black hole mass seems modest. When we con¬ 
sider that systematic errors may be contributing to the He II 
radial velocity width and that the HI radial velocity amplitude 
may be closer to the intrinsic amplitude, this only allows Mi 
to exceed 4 . 5 M 0 for i < 11° (probability 2%). Thus we con¬ 
clude that the modest radial velocity amplitude observed for 
IFGL 1018.6—5856 prefers a neutron star. However, these 
statistical statements do not exclude either option and a clean 
measurement of the companion radial velocity would help the 
situation greatly. 

The way forward on determining the mass of IFGL 
J1018.6—5856 doubtless starts with improved spectroscopy. 
Our data show that the photospheric lines are quite broad, 
and likely variable. Accordingly a future campaign on small- 
medium class telescopes should obtain long-slit, moderate 
resolution (R ~ 3000 — 6000) spectra covering 4000-5500A, 
with a collection of several template 06 V radial velocity 
comparisons. Queue/service observing is almost certainly re¬ 
quired in view of the large Pf,. Higher resolution echelle data 
from a large telescope can be of use, but will require higher 
S/N to allow deconvolution of the line shape and isolation of 
the systematic radial velocity. If the minimum mass exceeds 
3Mq, we can be confident of a black hole primary. But in 
the more likely event that lower masses are allowed, the incli¬ 
nation becomes critical and one will require very high preci¬ 
sion photometry to constrain the small ellipsoidal variations. 



Fig. 5.— The binary mass-mass plane. For the He Il-estimated radial veloc¬ 
ity amplitude, we see that a neutron star (NS) solution requires small compan¬ 
ion masses and/or large inclinations {i > 50°); roughly half of the allowed 
parameter space is consistent with a black hole (BH), although typical BH 
masses require very small inclinations. For the H I amplitude, a neutron star 
is highly preferred, requiring small inclinations to reach even IAMq. 
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or constraints from light curve fitting at other wavebands. 
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